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CHAPTER TWO 



ONE DIMENSIONAL WAVES IN SIMPLE SYSTEMS 
1 General solution to wave equation 



It is easy to verify by direct substitution that the most general solution of the one 
dimensional wave equation: 

^±= 2 ^± n n 

dt 2 C dx 2 1 ' } 

can be solved by 



where /(£) and g{£) are arbitrary functions of £. In the x, t (space, time) plane f(x — ct) 
is constant along the straight line x — ct = constant. Thus to the observer (x,t) who 
moves at the steady speed c along the positivwe x-axis., the function / is stationary. 
Thus to an observer moving from left to right at the speed c, the signal described 
initially by f(x) at t — remains unchanged in form as t increases, i.e., / propagates 
to the right at the speed c. Similarly g propagates to the left at the speed c. The lines 
x — ct =constant and x + ct = contant are called the characteristic curves (lines) along 
which signals propagate. Note that another way of writing (1.2) is 



= f(x 



ct) + g(x + ct) 



(1.2) 



<f>(x, t) = F{t - x/c) + G(t + x/c) 



(1.3) 



Let us illustrates an application. 



2 Branching of arteries 



References: Y C Fung : Biomechanics, Circulation. Springerl997 
M.J. Lighthill : Waves in Fluids, Cambridge 1978. 
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Recall the governing equations for pressure and velocity 



dt 2 " dx 2 
The two are related by the momentum equation 



d 2 u 9 d 2 u , n _ 

= ^ (2-2) 



The general solutions are : 



Since 



and 



p Tt = ~dx (2 ' 3) 



p = p+(x — ct) + p_(x + ct) (2.4) 
u = u + (x — ct) + u~(x + ct) (2.5) 



dp , , 

Wx =V ^ V - 



du- , , 

p-Qj- = ~pcu + + pcu__ 



where primes indicated ordinary differentiation with repect to the argument. Equation 
(2.3) can be satisfied if 

p + = pcu + , p_ = —pcu_ (2.6) 
Denote the discharge by Q = uA then 

ZQ ± = u±A = ±p± (2.7) 

where 

is the ratio of pressure to flux rate and is call the impedance. It is the property of the 
tube. 

Now we examine the effects of branching; Referening to figure x, the parent tubes 
branches into two charaterized by wave speeds C\ and C<i and impedaces Z\ and Z 2 . 
An incident wave approaching the junction will cause reflection in the same tube 

p = Pi(t-x/c)+p r (t + x/c) (2.9) 



2.2. WAVES DUE TO INITIAL DISTURBANCES 



4 



and transmitted waves in the branches are p\(t — x / c±) and p 2 (t — x/c 2 ). At the junction 
i = 0we expect the continuity of pressure and fluxes, hence 

Pi(t)+p r (t)=Pi(t)=p 2 (t) (2.10) 

^ = t + | ^ 
Define the reflection coefficient R to be the amplitude ratio of reflected wave to incident 
wave, then 

z ^ \z 1 ^ z 2 J 

Similarly the tranmission coefficients are 

T = Pi(i) kW I (213) 

Mt) P ,(t) i + + 

Note that both coefficients are constants depending only on the impedances. Hence the 
transmitted waves are similar in form to the incident waves except smaller by the factor 
T. The total wave on the incidence side is however very different. 

3 Waves in an infinite domain due to initial distur- 
bances 

Recall the governing equation for one-dimensional waves in a taut string 

d 2 u 2 d 2 u 

— -c— = 0, -oc <a ;<oc. (3.1) 
Let the initial transverse displacement and velocity be given along the entire string 

u(x,0) = f(x) (3.2) 

^(x,t) = g(x), (3.3) 

where f(x) and g(x) are non-zero only in the finite domain of x. At infinities x — > ±oo, u 
and du/dt are zero for any finite t. These conditions are best displayed in the space-time 
diagram as shown in Figure 1. 
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u -cu 

tt xx 



7 7 
u=f(x) u=g(x) 



Figure 1: Summary of the initial-boundary- value problem 

In (3.1.1) the highest time derivative is of the second order and initial data are 
prescribed for u and du/dt. Initial conditions that specify all derivatives of all orders 
less than the highest in the differential equation are called the Cauchy initial conditions. 

Recall that the general solution is 

u = 0(f) + ^(77) = 4>{x + ct) + - ct), (3.4) 

where and ip are so far arbitrary functions of the characteristic variables £ = x + ct 
and 77 = x — ct respectively. 

From the initial conditions we get 

u(x, 0) = <f>(x) + i]){x) — f(x) 
^(x,0) = c<P'(x)-c^(x)=g(x). (3.5) 

The last equation may be integrated with respect to x 



- f g(x')dx' + K, (3.6) 



1 Xo 

where K is an arbitrary constant. Now <fi and ip can be solved from (3.1.6) and (3.1.7) 
as functions of x, 

<j>(x) = 1 -lf(x) + K} + ±-J*g(x')dx' 



^x)= 1 -[f(x)-K}-±-J*g(x')dx', 
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t 



Range of 



(x,t) 




' c 



A 1 



x-ct 



,4 



Domain of 
dependence 



c \ 



x+ct 



X 



Figure 2: Domain of dependence and range of influence 



where K and xq are some constants. Replacing the arguments of by x + ct and of ip 
by x — ct and substituting the results in u, we get 



which is d'Alembert's solution to the homogeneous wave equation subject to general 
Cauchy initial conditions. 

To see the physical meaning, let us draw in the space-time diagram a triangle formed 
by two characteristic lines passing through the observer at x, t, as shown in Figure 2. 
The base of the triangle along the initial axis t = begins at x — ct and ends at x + ct. 
The solution (3.1.9) depends on the initial displacement at just the two corners x — ct 
and x + ct, and on the initial velocity only along the segment from x — ct to x + ct. 
Nothing outside the triangle matters. Therefore, to the observer at x, t, the domain 
of dependence is the base of the characteristic triangle formed by two characteristics 
passing through x, t. On the other hand, the data at any point x on the initial line t — 
must influence all observers in the wedge formed by two characteristics drawn from x, 
into the region of t > 0; this characteristic wedge is called the range of influence. 

Let us illustrate the physical effects of initial displacement and velocity separately. 





■x+ct 



(3.7) 
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Figure 3: Waves due to initial displacement 
Case (i): Initial displacement only: f(x) ^ and g(x) = 0. The solution is 

u{x,t) = ^f(x - ct) + -f(x + ct) 

and is shown for a simple f(x) in Figure 3 at successive time steps. Clearly, the initial 
disturbance is split into two equal waves propagating in opposite directions at the speed 
c. The outgoing waves preserve the initial profile, although their amplitudes are reduced 
by half. 

Case (ii): Initial velocity only: f(x) = 0, and g(x) ^ 0. Consider the simple example 
where 

g(x) = go when \x\ < b, and 
= when |x| > 0. 

Referring to Figure 4, we divide the x ~ t diagram into six regions by the characteristics 
with B and C lying on the x axis at x = — b and +b, respectively. The solution in various 
regions is: 

u = 



in the wedge ABE; 
in the strip EBIF; 



1 f x+ct g a 
u = — J g dx' = ^(x + ct + b) 



■y rx+ct 

u = — g Q dx = g t 

^ c Jx-ct 
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u, t 
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in the triangle BCI; 



in the wedge F/G; 



Figure 4: Waves due to initial velocity 



If,, So& 

« = 7T / 9odx = — 
2c i_ b c 



— / g dx' = —{b — x + ct) 
2c Jx-ct 2c 



in the strip GICH; and 



u = 



in the wedge HCD. The spatial variation of u is plotted for several instants in Figure 
4. Note that the wave fronts in both directions advance at the speed c. In contrast to 
Case (i), disturbance persists for all time in the region between the two fronts. 

4 Reflection from the fixed end of a string 

Let us use the d'Alembert solution to a problem in a half infinite domain x > 0. Consider 
a long and taut string stretched from x = to infinity. How do disturbances generated 
near the left end propagate as the result of initial displacement and velocity? 
At the left boundary x = must now add the condition 



u = 0, x = 0, * > 0. 



(4.1) 



In the space-time diagram let us draw two characteristics passing through x, t. For 
an observer in the region x > ct, the characteristic triangle does not intersect the time 
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axis because t is still too small. The observer does not feel the presence of the fixed end 
at x — 0, hence the solution (3.1.9) for an infinitely long string applies, 

fX + Ct 



1 1 r x + ct 

u = 7) [f( x + ct ) + f( x ~ ct )\ + 7T 0( r ) dr > x > ct - 

* ^ c Jx-ct 



(4.2) 



But for x < ct, this result is no longer valid. To ensure that the boundary condition 
is satisfied we employ the idea of mirror reflection. Consider a fictitious extension of 
the string to — oo < x < 0. If on the side x < the initial data are imposed such that 
f(x) = —f(—x), g(x) = —g(—x), then u(0,t) = is assured by symmetry. We now 
have initial conditions stated over the entire x axis 

u(x, 0) = F(x) and u t (x,0) — G(x) — oo < x < oo, 

where 



F(x) 



f(x) if x > 
-f(-x) if x < 



g(x) if x > 
— g(— x) if x < 0. 



G(x) = 

These conditions are summarized in Figure 5. Hence the solution for < x < ct is 



u = \ [F(x + ct) + F(x - ct)} + ^ ( f + / + ) G(x')dx' 

* ^ C \Jx-ct J0 J 

= \ [f(x + ct) - f(ct -x)] + y c (£ x + j^) 9( x ') dx ' 

1 1 r ct + x 

= - [f(x + ct) - f(ct -x)] + - g(x')dx'. (4.3) 

2 2c J ct _ x 

The domain of dependence is shown by the hatched segment on the x axis in Figure 6. 



5 Forced waves in an infinite domain 



Consider the inhomogeneous wave equation 

-QP = ° 2 Q^ + h &t) t > 0, \X\ < OO, 
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u=0 



u =c u 

tt xx 



u=-f(-x) 
u =-g(-x) 



u=f(x) 
u=g(x) 



Figure 5: Initial-boundary-value problem and the mirror reflection 



x+ct=x n +ct n , 

^ 1 t 



x-ct=x Q -ct Q 



'*(•*•» ^o) 



(x -ct ,0) o (ct -x ,0)(x +ct ,0) 



Figure 6: Reflection from a fixed end 
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where h(x, t) represents forcing. Because of linearity, we can treat the effects of initial 
data separately. Let us therefore focus attention only to the effects of persistent forcing 
and let the initial data be zero, 

du 



u(x,0) = 0, 
The boundary conditions are 



dt 



= 0, (5.1) 

t=o 



u — > 0, |x| — > oo. (5.2) 
Let the Fourier transform of any function f(x) and its inverse f(a) be defined by 

/oo 
f(x)e- iax dx (5.3) 
-oo 



1 



The transformed wave equation is now an ordinary differential equation for the transform 
of u(x, t), i.e., u(a, t), 

d U oo . _ 

— + c 2 a 2 u = h t>0 
at 2 

where h(a,t) denotes the transform of the forcing function. The initial conditions for u 
are: 

u(a, 0) = /(a) = 0, ^l°)=^(a) = 

Let us hide the parametric dependence on a for the time being. The general solution 
to the the inhomogeneous second-order ordinary differential equation is 

u = du^t) + C 2 u 2 (t) + J ^ [«i(r)fZ 2 (t) - u 2 (r)u 2 (t)} dr, (5.5) 
where U\ and u 2 are the homogeneous solutions 

«i = e- iact u 2 = e iact 

and W is the Wronskian 

W = u\u 2 — u 2 u x = 2iac = constant. 
The two initial conditions require that C\ = C 2 = 0, hence the Fourier transform is 



u 



r* H^r) iac{t _ T) _ e _ iac(t _ r) , dr (5 g) 

Jo 2*®c 1 J K J 
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The inverse transform is 

ft 1 POO 



U (x,t) = — [ dT- [ e -(x+c(t-r)) _ e ia(,-c(t-r))] rfa (5J) 



The integrand can be written as an integral: 

1 /** 1 Z" 00 /"X+c(t— 7") 



u = — 
2tt 



/ rfr— / dah(a,r) / d£e iaf 

Jo 2c y_ 00 } X - C (t- T ) 



By interchanging the order of integration 

^ /•< PX+c(t-T) ; 

' x—c{t~r) 

and using the definition of Fourier inversion, we get 



w = — 
2c 



/•t rx+c(t—T) i />oo 

/ rfr / / dae ia ^h(a,r) 

JO Jx-c(t-r) 27T 



u(x,t) = - / rfr / d£h(Z,T) (5.8) 

ZC JO Jx-c(t-r) 

The right-hand side is the integration of h over the characteristic triangle defined by 
the two characteristics passing through (x,t). in the x — t plane. Thus the observer is 
affected only by the forcing inside the characteristic triangle. 

For non-zero initial data u(x, 0) = f(x) and u t (x, 0) = g(x), we get by linear super- 
position the full solution of D'Alambert 

1 1 f*x~\~ct 

i(x,t) = -\f( x + ct) + f(x-ct)] + -J d£g(0 

i pt px+c(t-r) 

+ dr dCh(C,r), (5.9) 

Ac JO Jx-c(t-r) 

The domain of dependence is entirely within the characteristic triangle. 



u 



6 String embeded in an elastic surounding 

Reference : Graff : Waves in Elastic Solids 

If the lateral motion of the string is restrained by elastic springs along the entire 
length, the governing equation can be found from §1.1 by replacing the external pressure 
by the elastic restoring force — KV per unit length, 

d 2 V d 2 V 
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which can be written as 

1 d 2 V 

~c\d¥ ~~ 

where 

Co 

6.1 Monochromatic waves 

For any linearized wave problem, if the range of the spatial corrdinate x is (—00,00), 
and all coefficients are independent of x, t, then the first task is to examine the physics 
of sinusoidal wave train of the form: 

V(x, t) = \A\ cos(kx - iut -<j> A ) = U (Ae lkx ~ iu)t ) (6.4) 

where A = \A\e % ^ A is a complex number with magnitude \A\ and phase angle 4>a- After 
examining the physical meaning of this special type of waves, it is possible to use the 
principle of superposition to construct more general solutions. It is customary to omit 
the symbol 3?= "the real part of for the sake of brevity, i.e., 

V(x,t) = Ae lkx - iuJt (6.5) 

First of all a few definitions about sinusoidal waves in general. We shall call 

6(x,t) = kx - ut (6.6) 

the wave phase. Clearly the trigonometric function is periodic in phase with the period 
2tt. In the x, t plane, V has a constant value along a line of contant phase. In particular, 
6 = 2nn, (n = 0, 1, 2, • • • ) correspond to the wave crests where V = \A\ is the greatest. 
On the other hand, 6 = (2n + l)ir, (n = 0, 1, 2, • • • ) correspond to the wave troughs 
where V = — \A\ is the smallest. \A\ is half of the separation between adjacent crests 
and troughs and is called the wave amplitude; we also call A the complex amplitude. 
Clearly || represents the number of phase lines per unit distance, i.e., the density of 
phase lines, at a given instant; it is called the wavenwmber, 

wavenumber = k = — 
ox 



d 2 V K, 



dx 2 



-V 



(6.2) 



(6.3) 



(6.7) 
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On the other hand — 1| represensts the number of phase lines passing across a fixed x 
per unit time; it is called the wave frequency. 



wave frequency = uj = 



09 
"dt 



(6i 



To stay with a particular line of contant phase, say a crest, one must have 

dO = kdx — ujdt = 
namely one must move at the phase velocity, 



dx 


UJ 


c— ~r 
dt 


e=constant ^ 



(6.9) 



Now back to the string. Substituting (6.10) in (6.1) we get 

,2' 



, 2 K ur 
-k 2 - - + - 
J. ci 



V = 



or 



uj = c \ k 2 + 



K 
T 



The phase speed is 



UJ 

W\ 



K 
TP 



1/2 



> C n 



(6.10) 



(6.11) 



(6.12) 



See figure 6.1. Note that, due the stiffening by the lateral support, the phase speed 
is always greater than c G , and decreases monotonically with the wave number. Longer 
waves (small \k\) are faster, while shorter waves (larger \k\) are slower. As \k\ increases, 
c approaches the finite limit c Q for the shortest waves. In general a sinusoidal wave whose 
phase velocity depends on the wavelength, i.e., uj is a nonlinear function of k, is called 
a dispersive wave, and (6.11) or its equivalent (6.20) is called the dispersion relation. 

An interesting physical feature for dispersive waves in general can be found by su- 
perposing two trains of sighltly different frequencies and wave numbers: 



where 



V = A ^ e i ( k+x - UJ+t ) _|_ Ae^ k x ~ w *^ 



k ± = k±k', uj ± = uj(k ± ), k' < k. 



(6.13) 



(6.14) 
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Figure 7: Phase and group velocities of a string in elastic surrounding. Point of Sta- 
tionary phase k Q . 



Let us approximate by 



^ = Lj(k)±k ,d ^ + 0(k' 2 ) 
dk 



then 

y _ j^ e ikx-iut ^g-ifc' (x-(dw/dk)t) _|_ e ~ik' (x-(dw/dk)t)^j _ j^^. ^^kx-iut 

where 

.4 = 2cosA; / (a;- (dw/dk)t) (6.16) 

The factor exp(i/ca; — iut) is called the carrier wave and ^4(x,t) the envelope. Thus the 
result is a sinusoidal wave train with a slowly varying envelope which has a very long 
wavelength 2n/k' ^> 2n/k and moves at the group velocity 



group velocity = c g = — - (6.17) 

dk 



which is in general different from the phase velocity for dispersive waves. 

In our string problem, the group velocity is easily found from the dispersion relation 
(6.11) 



2 



c 9 = - C ° k = = ^ = —, hence °^ = %<l (6.18) 
9 ^k 2 + K/T c poo' c c 2 
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Thus the group velocity is always smaller than the phase velocity, and increases with 
the wavenumber from for the longest wave to the finite limit c G (equal to the phase 
velocity) for the shortest waves. 



When uj > uo c where 

ik 

T 

is called the cutoff frequency u c , k is real 



(6.19) 



fc = ±|fc|, with \k\=^—-—j (6.20) 

being the real and positive square root. V(x,t) is a propagating wave, with the plus 
(minus) sign corresponding to right- (left-) going wave. If u < u c , k = ±ik is imaginary, 
with k being the postive root: 

,2\ 1/2 



Then e tkx = e^ KX . For boundedness one must choose the minus (plus) sign for x > 
(x < 0). Oscillations are localized or evanecent; there is no wave radiation. 

As a simple application, (6.4) is the response in a semi-infinite string forced to 
oscillate at the left end x — 0, 

V(0,t) = ^t{Ae- iut }. (6.22) 

If uj > cu c , then 

V(x, t) = 3? (^e* 1 * 1 *-*"*) (6.23) 

where \k\ is defined by (6.20). The requirement that waves due to a local disturbance 
can only radiate outwards is called the radiation condition. More will be said on it in 
Appendix B. 

If uj < u c , we must require boundedness at infinity so that 

V(x,t) = $lAe- KX - iut , x>0 (6.24) 

At the cutoff frequency, k = 0, V is constant in x; the infinitely long string would 
oscillate in unison. This absurd result signfies the breakdown of the the linearized 
theory. 



2.6. STRING EMBEDDED IN AN ELASTIC SUROUNDING 



17 



6.2 Energy transport 

Along a unit length of the string the kinetic energy density is, 

The potential energy in any segment dx of the string is the work needed to deform it 
from static equilibrium. The part due to lengthening of the string against tension is 



fdV\ 2 T fdV\ 2 

TVls ,/,-) = /V 1 + dx -Tdx^- i^—j dx 



Adding the part against the springs, the total potential energy per unit length is 

We now calculate their time averages. If two time-harmonic functions are written in the 
complex form: 

a = 3? (Ae- iwt ) , 6 = 3? {Be~ iujt ) (6.27) 
the time average of their product is given by 

d>= — ! " abdt = \u{AB*) = -^(A*B) (6.28) 
2vr J t 2 2 

Using this recipe, the period-averaged energy densities are, 



KE = ^(-t^Ae^) 2 = ^lu 2 \A\ 2 



PE = ^$l (ikAe*) 2 + y (3L4e l9 ) 2 = {^k 2 + ^ \A\ 
where 6 = \k\x — cut is the wave phase. Since 



2 T fi2 K 



\rho V T 
we have 

P 



PE= r -uj 2 \A\ 2 
4 

Thus the total energy is equally partitioned between kinetic an potential energies, and 

P .21 A\2 



E = KE + PE = -uj \A\ (6.29) 
2 
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Now the averaged rate of energy influx across any station x can be calculated from 
the product of transverse component of the tension and the transverse velocity, 

This power input is transported from left to right by the wave. Hence the speed of 
energy transport is the group velocty. This result is quite general for many physical 
problems and is not limited to the springs. 



7 Dispersion from a localized initial disturbance 

The solution for monochromatic waves already shows that waves of different wavelengths 
move at different velocities. What then is the consequence of an initial disturbance? 
Since a general initial disturbance bounded in space can be represented by a Fourier 
integral which amounts to the sum of infinitely many sinusoids with a wide spectrum, 
we shall employ the tools of Fourier transform. 
In addition to the governing equation: 

1 d 2 V d 2 V KV 



c 2 dt 2 dx 2 T ' 
we add the initial (Cauchy) conditions 



oo < :r < oo, t > (7.1) 



r(*,0) =/(*), (7.2) 
Let us define the Fourier transform and its inverse by 

/OO I poo 

e- ikx f(x)dx, f(x) = — e lkx f(k)dk (7.3) 
-oo 27r J-oo 



The transforms of (7.1) is 



c 2 dt 2 T 
and of the intitial conditions are, 



1 d 2 V l2 - KV n ,„ A . 

-k 2 V-—, t>0 (7.4) 



F(M) = /(*), ^§^ = (7-5) 
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The solution for the tranform is 

V(k,t) = f(k) cos ut (7.6) 

where 



The Fourier inversion is 



u = Co^k* + ^ (7.7) 
1 f°° 

V( x ,t) = — dkf(k)e lkx cos cut (7.8) 
2vr J_ 00 

Any real function f(x) can be expressed as the sum of an even and an odd function 
of x. For simplicity let us assume that f(x) is even in x so that f(k) is real and even in 
k, then 

1 [°° 

V(x,t) — — / dkf(k) cos kx cosut 
n Jo 

which can be manipulated to 

i r°° 

V(x, t) = — & / dkf(k) (e ikx ~ iujt + e lkx+tuJt ) dk (7.9) 
2tt Jo 

The first term in the integand represents the right-going wave while the second, left- 
going. Each part correponds to a superposition of sinusoidal wave trains over the entire 
range of wave numbers, within the small range (k, k + dk) the amplitude is f(k)/2. The 
function f{k)/2 is called the Fourier amplitude spectrum. In general explicit evalua- 
tion of the Fourier integrals is not feasible. We shall therefore only seek approximate 
information. The method of stationary phase is particularly useful here. It aims at the 
asymptotic approximation of the integral 

J(t) = / F(k)e ltm dk (7.10) 

J a 

for large t. Let us first give a quick derivation of the mathematical result. Assume that 
F(k),(j)(k) are ordinary functions of k. If t is large, then as k increases along the path 
of integration both the real and imaginary parts of the exponential function 

cos(t<f>(k)) + % sm(t(f)(k)) 

oscillates rapidly between -1 to +1, unless there is a point of stationary phase k Q within 
(a,b) so that 

d ^ = 4>'{k ) = 0, a<k <b. (7.11) 
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Figure 8: Neighborhood of Stationary Phase 

Then important contribution to the Fourier integral comes only from the neighborhood 
of k . Near the point of stationary phase, we approximate the phase by 

( f ) (k)=^k )+ l -{k~k f^"{k ) + --- 

and the integral by 

/(*) » F{k )e lt ^ jf'exp (j(k - k ) 2 <P"(k )^J dk 

With an error of 0(l/t), we also replace the limits of the last integral by ±oo; the 

justification is omitted here. Now it is known that 

oo r~ 



1/2 / X X 

+ - , if k e (a,b), 

(7.12) 

where the sign is + (or — ) if <f)"(k ) is positive (or negative). It can be shown that if 
there is no stationary point in the range(a,b), then the integral I(t) is small 

I(t)=Opj), if k <£(a,b). (7.13) 



e 

V * 

It follows that 
I(t) = / F(k)e um dk w F{k )e lt ^ K)±w/A 
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Let us apply this result to the right-going wave 

1 f°° 

V + {x,t) = —UJ dkfitye^W-^dk (7.14) 



For a fixed x/t =constant, we have 



K 

(f)(k) — kx/t — u(k), with x/t = £, u = c \l k 2 + — (7-15) 

For an observer travelling at the speed x/t less than the maximum u'(k), i.e., c Q , there 
is a stationary point k Q at the root of 

x /t = u/(k ) = -^ = -^= (7.16) 

The root k increases from zero for long waves to the maximum c Q for very short waves. 
If the observer is faster than c , there is no stationary point. 
Since 

4>"{k ) = -uj"(k ) < 0, (7.17) 

the final result is that 



v+(x , t) . ^) KK) IJL_^-^ 



The transform f(k) depends on the specific profile of the initial disturbance. For exam- 
ple, if 

'M = (7 ' 19) 

which has an area S and characteristic width b, the Fourier transform is 

f(k) = Se- Wb (7.20) 

To a specific observer identified by the speed of travel x/t < c Q , the approximate 
result can be viewed as a simple harmonic wave train with wave number k Q , frequency 
ui{k ) and amplitude 

A- KK) I 2n (7 21) 

A -^rvj^lK) (7 - 21) 
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y -X* a^A^y- — — a^ .. 



\tL . 



Figure 9: Snapshots of wave field on the side x > 0. 

Since u/(fc) increases from with increasing k to the finite maximum T/pc = c Q , an 
observer faster than c G sees no waves. However any observer slower than c G is accom- 
panied by a train of progressive sinusoidal waves. The local wavelength k D is such that 
the group velocty matches the observer's speed. The faster the observer, the shorter 
the waves. If a snapshot is taken, then the shortest wave, whose phase velocity is the 
lowest are seen at the front which moves as fast as the crests of the shortest waves. The 
crests of longer waves advance much faster than the local envelope. Because / (k) is the 
greatest at k — 0, the envelope of the longest waves which stays near the source, is the 
biggest. The envelope of the shortest waves is lower in amplitude and spreads out with 
the wave front. The entire disturbance attenuates in time as t^ 1 ^ 2 . See figure 7 for an 
overview. 

Note also that very near the wave front, c g — * c Q ; the second derivative 4>"{k ) = 
—uj"{k ) vanishes. Hence the asymptotic formula breaks down. A better approximation 
is needed, but is omitted here. (See C. C. Mei, 1989, Applied Dynamics of Ocean Surface 
Waves). Finally we examine the propagation of wave energy in this transient problem. 
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Using (7.21) the local energy density is: 

F _ 1 „ .21 A |2 _ P^ KKf 

At any given t, the waves between two observers moving at slightly different speeds, 

Cg(ki) and c g (k 2 ), i.e., between two points x±/t — c g (k 1 ) and x 2 /t = c g {k 2 ) are essentially 

simple harmonic so that the total energy is 

x 2 r x 2 „,.,2 i fth \\2 



r%2 r 

I dxE= / 

J x-i J X 



4vrt w"(A; ) 
ui"{k )dk 



' X\ J X\ 

Since a; = u'(k )t for fixed t, we have 

dx 
T 

Now for x 2 > xi, k 2 > k±, it follows that 

rx 2 rk 2 2 

/ dxE = / dk Q — — (f(k )) 2 = constant 
J X1 Jki 47r 

Therefore the total energy between two observers moving at the local group velocity 
remains the same for all time. In other words, waves are transported by the local group 
velocity even in transient dispersion. 



8 Scattering of sinusoidal waves 

If along a long rod there are some localized inhomogeneities, an incoming train of sinu- 
soidal waves will be partly reflected and partly transmitted. The scattered signals tell 
us something about the scatterer. To determine the scattering properties for a known 
scatterer is called the scattering problem. To determine the scatterer from the scattering 
data is called the inverse scattering problem. We shall only study the former. 

Various mathematical techniques are needed for different cases: (i) Weak scatterers 
characterized by small amplitude relative to the wavelength, or slow variation within a 
wavelength, (ii) Strong scatterers if their dimensions are comparable to the wave length. 

8.1 Weak scattering 

Let the long rod have a slightly nonuniform cross section, 

S(x) = S (l + ea(x)), e < 1. (8.1) 
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where a(x) diminishes to zero at x ~ ±00. The wave equation reads 

-*!(«•>£) 

For sinusoidal waves 

u(x,t) = M [U{x)e~ iujt ] 
the spatial factor is governed by 

d ^S d f]+ P io 2 SU = (8.3) 



dx \ dx J 

Then the solution can be sought by the perturbation method 

U(x) = U (x) + eU^x) + e 2 U 2 {x) + ■■■ (8.4) 
Substituting this into the governing equation 

ES °ix~ { (1 + ea ^tx~ ^ o(x) + eUi ^ + e2{J ^ x) + " '1 } 

+pu 2 S (l + ea{x)) [U (x) + eU^x) + e 2 U 2 (x) + •••]= (8.5) 

and equating the coefficients of like powers of e to zero, we get from the terms of order 
0(e°): 

ES o d ^ + puj 2 S U = (8.6) 
The solution is simply the incident wave 



U = Ae ikx , with k = u^p/E (8.7) 



At the order 0(e) we get 



ES °~Q^2 + P UJ S ° Ul + a ( X ) \ ES °~Q^2 + PU S ° U ° ) °dx~~dx~ = 

or, 

^ + ^1 = -^ (8-8) 
Thus the scattered wave is at most of the order e and is governed by an inhomogeneous 
equation. Let us define the fundamental solution (Green's function) by 

d2G \ x ^ + k 2 G(x, x') = 5(x - x') (8.9) 
dx 2 
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G ~ outgoing waves at ± oo 
It will be shown in Appendix A that 

^ik\x—x'\ 



G(x,x') 



2ik 



no) 



Ul) 



By straightforward differentiation, it can be verified that 

da(x') dUo(x') 



ax dx 2ik l_ r ^ dx dx 



±A e ik\x- x >\ dx , (gJ2) 



Since U (x') = e lkx> , we have 
A Z" 00 da(x') 



2 ./_«, dx' 



Jkx' e ik\x-x'\^ x l 



A 



Akx 



da 
dx' 



dx + e" 



' — oo 

A ( 



f 

J X 



00 rln 

dx' 



(8.13) 



2 



e lkx a(x) + e 



-ikx 



^_ p 1ikx' j„/ 



ifcx i —ikx 



a(x)e lkx + e 



A 

dx 



e"""" dx 1 

- (a(x')e 2ikx '^ dx' - a2ike 2ikx ' 



dx' 



= ikAe~ ikx / a(x')e 2ikx ' dx' 

J X 

Far to the right we get 



(8.14) 



U! = 0, x^oo. (8.15) 

Thus on the transmission side modification of the incident waves can at most be of the 
order 0(e 2 ). Far to the left x ~ — oo, 



U**eU 1 = Ae~ ikx (iek j a(x')e 2ikx ' dx'^j , x 



-oo. 



Let us define the reflection coefficient by 

/oo 
a(x')e 2ikx ' 'dx' , 
-oo 

For instance if a is a hump of unit total dimensionless area and length b, 

a(x) = 



tt(x 2 + 6 2 ) 



then 



R = iekbe 



-2kb 



(8.16) 



(8.17) 



(8.18) 



(8.19) 
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Figure 10: Reflection coefficient given by ( 8.19). 



Hence reflection is small if kb is large, i.e., long and gentle obstacles are not effective 
in reflecting short waves . For very long waves or a short obstacle, reflection is also 
vanishingly small. The maximum reflection coefficient is R = iea a /2 when kb = 1/2. 
Note that far to the right, U\(x) — > as x — > oo. The transmission vwave is 



U^U + eU x = Ae ik ° x {\ + 0(e 2 )). 



(8.20) 



The tranmission coefficient is T = 1 + 0(e 2 ). This is consistent with the general law of 
energy conservation |i?| 2 + |T| 2 = 1 to be proven later. 

Homework: Check if the governing equation at 0(e 2 ) for U 2 is: 



ES 



d 2 U 2 , 2 



dx 2 



+ pufSoU* = ES 



dadU dadUi 



dx dx dx dx 



.21) 



Carry out the solution to find the transmission coefficient T from the amplitude of 
Uq + eUi + e 2 f/ 2 at x ~ oo, and show that energy is conserved to the order 0(e 2 ). 
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8.2 Strong scattering of long water-waves by a shelf 

Consider an ocean bottom with a step-wise variation of depth. 

hi, x < —a; 

h = { //•>. -a < x < a; (8.22) 
h s — hi, x > a 

How does the step interupt the propagation of an incident wave of unit amplitude 
arriving from x ~ — oo? 

In each zone of contant depth (i = 1, 2, 3), the shallow water equations read: 

a "- +9 t-0 (8.24) 



For sinusoidal waves 



we get 



dt dx 



Ci = ?7ie , u i = U i e- vM (8.25) 



-iurji + h~^ = (8.26) 



dx 

-iuUi + g^^O (8.27) 



or 



r2 

ar]i + kfa = 0, where h t = (8.28) 



dx 2 1 11 ' 1 Jgh, 

At a junction, the pressure and the flux must be equal, hence 

771 = 772, and hi^- = h 2 C ^, x = -a; (8.29) 
ax ax 

drj 2 dr] 3 

7/2=7/3, and h 2 — = hi — , x = a. (8.30) 
The forms of the solutions in each zone of constant depth are: 

771 = e ikl{x+a) + Re- ikl(x+a \ x < -a; (8.31) 

m = Ae ik2X + Be- ik2X , - a < x < a (8.32) 

773 = Te lkl{x ~ a \ x>a (8.33) 
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The reflection and transmission coefficients R and T as well as A and B are yet unknown. 
Applying the matching conditions at the left junction, we get two relations 

1 + R = Ae~ ik2a + Be ik2a 

hhxil -R) = k 2 h 2 {Ae- lk2a - Be ik2a ). 
Similarly the matching conditions at x = a gives 

Ae ik 2 a + Be ~ik 2 a = j, 

k 2 h 2 {Ae ik2a - Be~ ik2a ) = k x h x T. 
These four equations can be solved to give 

4s 



R 



'~y _|_ 5^2g— 2ik 2 a M g\2^2ik 2 a 

— (\ — s 2 )(e~ 2ik2a — e 2ifcaa ) 

'"Y _|_ £^2^—2ik 2 a g^2^2ik 2 a 

A = ^ e - ik2a (1 + s) 

B = -e ik2a (1 - s) 
2 v ; 

where 

k\h\ hi c\ 



k 2 h 2 V h 2 c 2 
The energy associated with the tranmsitted and reflected waves are : 

T 



\R\ 



4s 2 + (1 - s 2 ) 2 sin 2 2fc 2 a 
(1 - s 2 )sin 2 2A;2a 



4s 2 + (1 - s 2 ) 2 sin 2 2fc 2 a 
It is evident that \R\ 2 + \T\ 2 = 1. 

Over the shelf the free surface is given by 

2s [(1 + s )e ik2{x ~ a) + (1 - s)e _ifca ^-°)] 
V2 = (1 + s) 2 e- ik * a - (1 - s) 2 e ik * a ^ 8 ' 45 ^ 

Recalling the time factor e~' tuJt , we see that the free surface over the shelf consists of two 
wave trains advancing in oppopsite directions. Therefore along the shelf the two waves 
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can interfere each other constructively, with the crests of one coinciding with the crests 
of the other at the same moment. At other places the interference is destructive, with 
the crests of one wave train coinciding with the troughs of the other. The envelope of 
energy on the shelf is given by 

2 4s 2 [cos 2 k 2 (x — a) + s 2 sin 2 k 2 (x — a)] 
M = 4 S 2 + (l- s )2 sin 2 2 A; 2 a (8 ' 46) 

At the downwave edge of the shelf, x = a, the envelope is 

2 

|r?l = 4s 2 + (l-s) 2 sin 2 2A; 2 a (8 ' 4?) 

Note that the reflection and transmission coefficients are oscillatory in k 2 a. In par- 
ticulatr for 2k 2 a = nir,n = 1,2,3..., that is 4a/A = n, \R\ = and |T| = I. The 
shelf is transparent to the incident waves, correspoinidng to the most constructive inter- 
ference and the strongest transmission Mininum transmision and maximum reflection 
occur when 2k 2 a = (n — l/2)n, or 4a/A = n — 1/2, when the interference is the most 
destructive. The corresponding transimssion and reflection coefficients are 



min|T| 2 = 7-^^, max|i?| 2 = ^ (8.48) 



As 2 2 _ (1 -g 2 ) 2 

[iT^Y' max| 1 ~ fi + s 2 ^ 2 ' 



See Figure 11 

The features of interference can be explained physically. The incident wave train 
consists of periodic crests and troughs. When one of the crests first strikes the left 
edge at x — a, part of the it is transmitted onto the shelf and part is reflected towards 
x ~ —00. After reaching the right edge at x — a, the tranmitted crest has a part 
reflected to the left and re-reflected by the edge x = —a to the right again. When the 
re-reflected crest arrives at the right edge the second time, its total travel distance is 
4a. If 4a is an integral multple of the wave length A 2 , the crest is in phase with all 
the other crests entering the shelf either before (the third, fourth, fifth, ... time) or 
after. Thus all the crests reinforce one another at the right edge. This is constructive 
interference, leading to the strongest tranmission to the right x ~ 00. On the other 
hand if 2k 2 a = (n — l/2)ir or 4a/A = n — 1/2, some crests will be in opposite phase to 
some other crests, leading to the most destructive interference at the right edge, and 
smallest transmission. 
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Figure 11: Transmission and reflection coefficients for a rectangular shelf. From Mei, 
C.C., Applied Dynamics of Ocean Surface Waves. 

From the results in the following exercise, the reflected wave suffers a reversal of 
phase at the discontinuity (a crest turns into a trough and vice versa). Hence at the 
reflection edge a wave crest tranmitted onto the shelf returns from trough. 
It is then transmitted to the left of x = —a as trouigh. On the other hand, a crest that 
is reflected at x — —a appears as a crest. Therefore the interference is destructive. The 
resulting R, which is the sum of multiple reflections, is small. 

Homework: Scattering by a depth discontinuity: Consider the scattering of long 
water waves by a depth discontinuity at x — 0, i.e., h = hi,—oo < x < and h = 
h 2 , < x < oo. Show that the tranmssion and reflection coefficients are: 

T = 2 fi. llMlg (8 49) 

l + (/i 2 //n)V2' n i + ( h2 / hl y/2 V-™) 

Thus the reflected wave suffers a reversal of phase when h 2 /hi > 1, i.e, approaching the 
discontinuity from the shallow side. 
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9 General identities in scattering of monochromatic 
waves by arbitrary inhomogeneities 

Scattering due to inhomogeneities, caused by nonuniformities either in geometry or in 
material properties, requires the solution of ordinary differential equations with variable 
coefficients. In general one must resort to numerical means. For the sake of checking 
numerical accuracy and for gaining physical insight, identities which must be true are 
useful. They are often deduced by arguments typical in the derivation or the use of 
Green's theorem. 

We illustrate these identities by the example of an infinitely long rod with variable 
cross section S(x), 

dx ^dx) = P ~EW (9,1) 
E and p are taken to be constants. S(x) is nonuniform only in a finite neighborhood 
around x — 0. Elsewhere S = S Q =constant. Consider monochromatic waves 

u(x,t) = U{x)e~ iujt (9.2) 

so that U (x) is governed by the ordinary differential equation 

(SU')' + k 2 SU = 0, k 2 = ^f (9.3) 

E 

The boundary conditions depend on the source of the incident waves. 

Let us consider two problems: one with the incident wave from x ~ — oo; the other 
with the incident wave from x ~ oo. 

If the incident waves arrive from x ~ — oo, the asymptotic form of the solution must 

be 

U^x) ~ A„e ikx + B-e~ ikx = A. (e ikx + R ie - ikx ) , kx ~ -oo (9.4) 

and 

U ± {x) ~ A + e lkx = A_T ie ikx kx ~ oo (9.5) 

where 7\ and R\ are the transmission and reflection coefficients which are a part of the 
unknown solution of the left-incidence problem. Let us define the Jost function by 
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then 



/l(l) ~^ ib + | e_ih hx4H-\ (9.7) 



and 

/i(ar) ~ A;x > 1 (9.8) 

On the other hand, if the incident waves arrive from x ~ oo, the asymptotic form of 
the solution must be 

U 2 {x) ~ 6_e- fa = 6+T 2 e- ifex , A;x ~ -oo (9.9) 

and 

U 2 (x) ~ a+e^ + fr+e"*** = 6+ (i^e*** + e" ifca; ) , kx ~ oo (9.10) 

where T 2 and i?2 are the unknown transmission and reflection coefficients of the right- 
incidence problem. Similarly we define another Jost function 

h = T~^ (9.11) 

so that it behaves as 

f 2 (x)~e~ ikx , kx^-oo (9.12) 

and 

f2(x)~^e- lkx + I ^e* kx , kx~oc (9.13) 

J-2 J-2 

Since both f\ and f 2 satisfy (9.3), we have 

(Sf[)' + k 2 Sf 1 =0, and (Sf£ + k 2 Sf 2 = 0, (9.14) 

Multiplying the first by f 2 and the second by fx, then taking the difference, we get after 
partial integration, 

^ (su\f 2 - / 2 /0) = o, 

implying the following Wronskian identity 

W(f 1 J 2 ) = (f i n-f 2 f[) = ^- ) (9.15) 
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The constant C can be related to the asymptotic limits at k\x\ 3> 1 where S — > So- Far 
on the incident side, the Wronskian is 

(/1/2 ~~ hf[)x — 00 

1 ikx 1 ^1 -ib| / -u\„-ikx -ikx ( ^ Akx ikR\ ih x 



e lkx + _L e - tteB (-ik)e~ lkx - e~ lkx —e tkx -e 

2ik 

Far on the transmission side, the Wronskian is 

(/1/2 ~ f2f'i)x~oo 

= jk* ( Z^ e -ikx + ikR2 ikx\ _ ( 1 -ikx + R2 ikx\ , ik) ikx 

\T 2 T 2 J \T 2 T 2 ) y ' 

2ik 

Since the two Wronskians are equal, we conclude that 

T x = T 2 (9.16) 

Thus, the complex transmission coefficients are equal for incident waves coming from 
either side, no matter how asymmetric the geometry may be in the near field! 

Next we consider the left-incidence problem and omitting the subscript " 1" . Multi- 
plying the first of (9.3) by the complex conjugate of U, i.e., U*, and taking the difference 
of the result with its complex conjugate, we get 

US(U*J - U*(SU')' = ^- (S(UU*' - U*U')) = (9.17) 

hence 

S(UU*' - U*U') = constant (9.18) 
By using the asymptotic expression (9.4) on the incidence side, we get 

[S{UU*' -SUW)]^^ 

= S (A_e tkx + B_e~ ikx ) (-ikAW tkx + ikB*_e tkx ) 
-S [A*_e~ ikx + B*_e tkx ) (ikA„e lkx - ikB„e~ lkx ) 
= -2ik\A\ 2 + 2ik\B\ 2 
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Similarly by using the asymptotic expression (9.5) on the transmission side, we get 

[s{uu*' -u*u')] x ^ 

= S (A + e lkx ) {(-ik)A\e~ lkx ) - S Q (A* +e - lkx ) {(ik)A + e ikx ) 
= *2ikS 1 A_^_ | 

Equating the two limits we conclude that 

\A_\ 2 = \B_\ 2 + \A + \ 2 (9.19) 

or 

\R 1 \ 2 + |7\| 2 = 1. (9.20) 

This is merely a statement of conservation flux : the flux rate of total scattered energy 
(reflected and transmitted) is equal to the flux rate of the incident wave energy; the 
speed of energy tranport being the same on both sides of the scatterer. Clearly the 
same energy conservation must hold for the right-incidence problem. 

\R 2 \ 2 + \T 2 \ 2 = 1 (9.21) 

If the scattering problem is solved numerically by, say, finite elements, it is necesssary 
that the computed scattering coefficients R and T satisfy the identities (9.16) (9.20) and 
(9.21). 

Homework: Extend the identities (9.16) and (9.20) when S(x) — > Si as kx — > — oo 
and S(x) — > S2 as kx — > 00, where Si and S2 are two different constants. 

10 Refraction in a slowly varying medium 

Let us consider time-harmonic long waves in shallow water waves. The governing equa- 
tion is 

Consider a sea depth which varies slowly within a wavelength, i.e., 
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Earlier analysis suggests that reflection is negligibly small. Thus the solution is expected 
to be a locally progressive wave with both the wavenunmber and amplitude varying much 
more slowly than the wave phase in x . Hence we try the solution 

r] = A(x)e i6 ^ (10.3) 

where 9{x) — out is the phase function and 

d9 

k(x) = 7- (10.4) 
ax 

is the local wave number and A is the complex amplitude. Note the spatial rate of 
variation of the phase function, i.e., the wave number is in gneral not small, therefore 9 
itself is not a slowly varying function of x. Let us calculate the first derivative: 

and assume 

dA 

-fr = O(kL)- 1 < 1 

In fact we shall assume each derivative of h, A or k is \i times smaller than kh, kA or 
k 2 . Futhurmore, 



d f drj\ uj 



,2 



... . _ . . , <l.\\ d ,<l.\\ d(khA) u?A ... 

, ih-r)+—v= ik ikh + h— + — h— + i v . ' + e w = 

dx \ dx J g \ dx J dx \ dx j dx g 

Note that the complex amplitude A can be written as 

A = \A{x)\e tdA(x) (10.5) 

The phase of A, 9 A) is a slowly varying function of x. It can be considered as a part of 
the wave phase, although its spatial gradient is much smaller than k. 
Now let's expand 

A = A + A 1 + A 2 + --- (10.6) 
with Ai/A = 0(/x), A 2 /A = 0(/x 2 ), • • • . From 0(/jP) the dispersion relation follows: 

u 2 = ghk 2 , or k = (10.7) 

Vgh 
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Thus the local wave number and the local depth are related to frequency according 
to the well known dispersion relation for constant depth. As the depth decreases, the 
wavenumber increases. Hence the local phase velocity 



also decreases. 

From O(n) we get, 

or, after multiplication by A , 

which implies 

hence 



c=^ = ^Jh. (10.8) 



ax ax 



d (kh\A \ 2 ) = (10.9) 



dx 



or, 



kh\A \ 2 = C 2 = constant 



\fgh\A \ 2 = constant = a/^ooI^ooI 2 (10.10) 



Since in shallow water the group velocity equals the phase velocity, the above result 
means that the rate of energy flux is the same for all x and is consistent with the 
original assumption of unidirectional propagation. Furthermore, the local amplitude 
increases with depth as 

A Q (x) /// x n 



Aqq \ h 

This result is called Green's law. 

In summary, the leading order solution is 

V4 /U \ V4 



(10.11) 



C = Aoo (^j e ie -^ = A^ (^j exp (i J X k(x')dx' - iu?J (10.12) 

Homework: Study the refraction of a train of obliquely incident waves towards a 
sloping beach. The sea depth varies slowly only in x, h — h(x). 
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A Solution for the Green function 



Let us divide the domain into two parts, corresponding to two sides of the concentrated 
forcing. On the right side x > xo, we define G(x,xo) = G+(x, xq) which must be 
governed by : 



d 2 G 
dx 2 



+ k^G + = 0, x > x 



(A.l) 



To satisfy the radiation condition we take the solution to beTo satisfy the radiation 
condition we take 

G + = Ce ikoX (A.2) 
On the left side we define G(x,x ) = G-(x,x ), which must be governed by 

d 2 G- 



dx 2 



+ k G- = 0, x < xq 



(A.3) 



To satisfy the radiation condition we take 

G + = De~ ik ° x 

At x = Xo we require continuity: 

G + (x ,x ) = G-{x),X)) 

Hence 



(A.4) 



Ce 



ik XQ 



De 



-ikoxo 



(A.5) 



(A.6) 



In addition we get another matching condition at x = Xo by integrating eq(8.8) across 
the concentrated forcing 

rxo+8 r (JZq ~\ rxo+5 

/ < — — - + k^G > dx = / S(x — x )dx = 1. 

Jxo-S {.MX J Jx {) -8 



The left-hand side can be integrated to give the second matching codition, 

= 1 



'dG+ 




~dGJ 


dx 


XQ 


dx 



implying that 



ikoCe ikoxo + ik De- iki,x ° = 1 



(A.7) 



(A.8) 
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Thus 



The final solution is 



c ° = w Do= w (A ' 9) 



e iko(x-xo) e ~ik (x-x ) 

G + = —- , x-x >0; G- = — — , x-x <0. 

or more compactly, 

„iko\x— xo\ 

G(x,x )= (A.10) 
More on one-dimensional waves will be discussed later. 



B On the radiation condition 
B.l Forms of the radiation condition 

In the theory of diffraction and/or radiation of time-harmonic waves, an important 
boundary condition is the so-called radiation condition due to A. Sommerefeld : 

The disturbance due to either diffraction by a finite scatterer, or radiation from a finite 
region, can only behave as outgoing waves at infinity. 
Let us first write the unknown as 

$(x,t) = 0(x)e^ (B.I) 

In a scattering (diffraction) problem we let 

= 0/ + 0' (B.2) 

where 0/ represents the incident wave (e.g., 0/ = A e lkx ), and 0' denotes the distur- 
bance (scattered wave). In a radiation problem, 0' represents the radiated wave due to 
motion in the finite region. The radiation condition can be stated mathematically in 
two equivalent forms, as follows. 

1. One-dimensional propagation: 

<j>'(x) -> A B±e ±tkx , as kx -> ±oo; (B.3) 

or, 

lim \^r- =F ik<p) -> (B.4) 

kx^±oo \ OX I 
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2. Two-dimensional propagation: 

pikr 

4>'(r,9,z) -> A B(6,z)—=, as kr -> oo; (B.5) 

v fcr 

or, 

lim ( ^ - ] -> (B.6) 

fcr^oo \ <7r / 

3. Three-dimensional propagation: 

0'(r,M) -> 4>L(M)-j— , as A;r^±oo; (B.7) 

or, 



lim fcr - ik<j>' -> (B.8) 
fcr->oo y ar / 

Note that the radiation condition is meant only for waves that is strictly sinusoidal in 
time from ancient past. For transient disturbances, the condition should be changed to 

$'(x, — > 0, as |x| — > oo, t < oo. (B.9) 

Physically, if the source of the wave disturbance is in a finite region, it is not possible 
for wave energy to flow inward from infinity. Mathematically, a sinuoidal wave is the 
limit of a transient wave where the local and sinusouidal forcing begins at some time 
in the past, say t — 0. For all finite t, the boundary condition (B.9) holds. In order 
to derive the behavior at very large time (the quasi steady state), one can solve an 
initial-bounday value problem under the initial conditions, 

0',^ = O, t = 0, Vx (B.10) 

and the additional boundary condition in the source region, e.g., 

^ = n-Ve- !Wt , t>0. (B.ll) 
on 

and then take the limit of large t. The explicit solution for the IBVP is not trivial. We 
shall demonstrate for a one-dimensional example where the local energy supply grows 
slowly in comparison with the wave period. 



2.10. REFRACTION IN A SLOWLY VARYING MEDIUM 



40 



B.2 Transient motion of an elastically supported string 

Consider a semi-infinte string, 

pV tt - TV XX + KV = 0, , x>0, t>0 (B.12) 

subject to the boundary condition, 

V(0,t) = A (et)e- iojt , t > (B.13) 

where e C 1 and A (et) increases slowly from A (— oo) = to A (oo) = A =constant 
with the time range ujt\ = 0(1). Thus at the boun deary there is sinusoidal forcing with 
a slowly increasing amplitude. At infinity the boundary condition is 

V(oo,t)->0, t<oo. (B.14) 

Let us use the perturbation method of mulitple scales. In vieew of (B.13), we in- 
troduce two time coordinates t and t\ = et. Suggested by the elementary example of 
wave groups, we expect that the envelope to be slowly varying in space as well. Let 
us therefore also introduce two space coordinates x, x\ = ex±. We now pretend that 
V = V(x,Xi,t,ti) so that the original time derivatives become, according to the chain 
rule, 

d dV dV 

m^lH +e W 1 > (R15) 

^ ^ + e|U? + ^^ + 2 ^ + ^. (B.16) 



dt 2 \dt Oh J \dt dtj dt 2 dtdh dt\ 
Likewise the space derivatives become 

d ^dV dV 

dx dx dx 1 ' 

d 2 V (d d\fdV dV\ d 2 V d 2 V 2 d 2 V . 



dx 2 \dx dx 1 J \dx dx 1 J dx 2 dtdx 1 dx\ 
We next substitute these into (B.12) 

V = Vo(x,Xi,t,ti) +eV 1 (x,x 1 ;t,t 1 ) + e 2 V(x, x^t,^) + ■■■ (B.19) 

and separate different orders. At order O(e ) we have the following perturbation equa- 
tion : 
0(e°): 

P^-1& + KV a = (B,0, 
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At order 0(e), we get: 

Let us take the leading-order solution to be a sinusoidal progressive wave 

V = A{x 1 ,t 1 )e i( * x - ut) (B.22) 

where 

* = ±V*f-f (R23) 

Let the square root be real so that there is propagation. However, we do not know the 
sign of k, hence the direction of propagation. 
At the order 0(e) Eqn. (B.21) can be written 

+ W = -» f k Tl± + ^ . (B.24) 



<9£ 2 <9x 2 \ 9xi 

The forcing term on the right is a homogeneous solution of the left, hence unbounded 
resonance of V\ must occur unless 

dA OA 

*W""toi= ' (R25) 

or 

This is the condition of secularity (solvability) and is a first order hyperbolic equation. 
The general solution is 

^=4'-^) (R27) 

Now we apply the boundary condition (B.13) at x\ = 0, 

A{U) = Aoih) (B.28) 

Hence the envelope is 

A ^ = A '{ h -W^) < R29 > 

For this transient problem we expect that A(oo,ti) approaches zero when x — ■> oo 
for finite ti, according to (B.14). Since for large x±, A(x 1 ,ti) — > A (—xi/Tk/puj), we 
must take the positive sign for k in view of the given behaviour of A . Thus 

k > (B.30) 
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implying that the solution is a rightward propagating wave. Thus the radiation condition 
is derived. 

Note further that 

Tk 

— =c 9 (B.31) 

puj 

is the group velocity in the string, thus 

A(x 1 ,t 1 ) = A(h -x 1 /c g ) (B.32) 
or, A is a slowly varying function of x — c g t. Thus result is true for all dispersive waves. 



